Experimental demonstration of single-site addressability 
in a two-dimensional optical lattice 
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We demonstrate single site addressability in a two-dimensional optical lattice with 600 nm lattice 
spacing. After loading a Bose-Einstein condensate in the lattice potential we use a focused electron 
beam to remove atoms from selected sites. The patterned structure is subsequently imaged by means 
of scanning electron microscopy. This technique allows to create arbitrary patterns of mesoscopic 
atomic ensembles. We find that the patterns are remarkably stable against tunneling diffusion. Such 
micro-engineered quantum gases are a versatile resource for applications in quantum simulation, 
quantum optics and quantum information processing with neutral atoms. 



PACS numbers: 03.75.Hh, 37.10.Jk, 62.23.St 

Ultracold atomic gases provide a versatile experimen- 
tal platform to address fundamental quantum mechanical 
phenomena. On the one hand, they are a subject of re- 
search on their own [TJ [2]. On the other hand, they are 
a resource for applications in quantum optics [3], quan- 
tum simulation [4] or quantum information processing 
[5j E]. In order to use neutral atoms for such purposes, 
micro-structured quantum gases with suitable addressing 
schemes are essential. Optical tweezers [HE] [9], optical 
conveyor belts [10] and magnetic microtraps [TTJ [12] are 
currently investigated in this context. Ultracold quantum 
gases in optical lattices are an alternative as they can 
combine perfectly periodic potential structures in one, 
two and three dimensions with single site occupancy [13] . 
However, a lattice spacing below 1 fim imposes stringent 
requirements on the spatial resolution of the addressing 
scheme. Several approaches are currently pursued in or- 
der to tackle this difficulty. Advanced optical addressing 
schemes are a natural choice [14j [15] . Coupling atoms to 
single ions [16] or a small solid probe are alternatives to 
reach this goal. Lattices with larger spacing [17] are also 
promising as they relax the requirements for the nominal 
resolution of the addressing scheme. 

Here, we demonstrate single site addressability of a 
quantum gas in a 2D optical lattice with 600 nm lattice 
spacing. To this end, we employ a recently developed 
imaging and manipulation technique based on scanning 
electron microscopy [18] . Loading a Bose-Einstein con- 
densate in the lattice potential, we selectively remove 
atoms from individual sites by means of the dissipative 
interaction with a focused electron beam. In this way, ar- 
bitrary patterns of occupied lattice sites can be produced. 
We study the tunneling processes in the patterned struc- 
ture and the depletion dynamics during the preparation 
stage. 

In our experiment we use a scanning electron micro- 
scope which is implemented in a standard apparatus for 
the production of ultracold quantum gases [I9j [20] . Using 
all-optical techniques we produce a Bose-Einstein con- 




FIG. 1: Electron microscope image of a Bose-Einstein con- 
densate in a 2D optical lattice with 600 nm lattice spacing 
(sum obtained from 260 individual experimental realizations) . 
Each site has a tube-like shape with an extension of 6 /mi per- 
pendicular to the plane of projection. The central lattice sites 
contain about 80 atoms. 



densate of 10 5 rubidium atoms in a single beam optical 
dipole trap. The 2D optical lattice is generated by two 
independent laser systems. For each axis, a laser beam 
is split into two beams which interfere under an angle of 
90° at the position of the dipole trap. In this geometry, 
the laser wavelength of 850 nm results in a lattice spacing 
of 600 nm in both directions. The condensate is adiabat- 
ically transferred into the ground state of the combined 
optical potential ramping up the laser power of the lat- 
tice beams in 200 ms. We can achieve a maximum lattice 
depth of 18 recoil energies, and the filling factor in the 
central part of the gas is 80 atoms per site. Each site 
has a tube-like shape with an extension of 6 /am along 
the direction of the electron beam. During the imaging 
procedure the electron beam is scanned across the cloud. 
Electron impact ionization produces ions which are de- 
tected by a channeltron. The overall detection efficiency 
is about 15 %. Figure[T] shows a scanning electron micro- 
scope image of a Bose-Einstein condensate in the lattice 
potential. The image was obtained by summing over 260 
individual experimental realizations. 

The addressing of individual sites in the optical lattice 
works as follows: the electron beam is pointed at selected 
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FIG. 2: Patterning a Bose-Einstein condensate in a 2D optical 
lattice with a spacing of 600 nm. Every emptied site was illu- 
minated with the electron beam (7nA beam current, 100 nm 
FWHM beam diameter) for 3 ms (a,b), 2 ms (c,d), and 1,5 ms 
(e), respectively. The imaging time was 45 ms. Between 150 
and 250 images from individual experimental realizations have 
been summed for each pattern. 



sites for a dwell time of 1-3 ms per site in order to remove 
the atoms. Immediately after the preparation, the imag- 
ing procedure is started. Four elementary examples of 
this patterning technique are presented in Fig. [2] A sin- 
gle defect in the lattice structure is shown in Fig.|2|a). 
The structure resembles a Schottky defect in a solid and 
is an ideal starting point to study the tunnelling dynam- 
ics close to a defect. The opposite situation corresponds 
to an isolated lattice site and is shown in Fig.[2jb). Such 
a mesoscopic ensemble provides, for instance, the possi- 
bility to study the transition from few-body systems to 
the thermodynamic limit. It can also act as a paradigm 
for Rydberg blockade studies as the spatial extension of 
the ensemble is very small. In this context, quantum op- 
tical applications such as single atom and single photon 
sources or the creation of GHZ-like states have been pro- 



posed 
in Fig. 



2Tj . A chain and a ring of lattice sites are shown 
2^c) and [2^d) in order to illustrate the large va- 



riety of achievable geometries. Obviously, our approach 



allows for any arbitrary pattern that fits to the underly- 
ing quadratic lattice geometry (Fig.^e)). 

Avoiding crosstalk to neighboring sites during the ad- 
dressing procedure imposes stringent requirements on the 
spatial resolution. This is in contrast to a pure imaging 
process, for which it has been shown that neighboring 
sites can be distinguished even with a resolution limit 
well above the lattice spacing [10]. Hence, imaging lattice 
sites and manipulating lattice sites set in general differ- 
ent requirements on the spatial resolution. As in our case 
the diameter of the electron beam (100 nm FWHM) is 
much smaller than the lattice spacing (600 nm), the sites 
can be addressed without affecting neighboring sites (see 
Fig.[2|a)). The presented scheme is especially attractive 
in combination with single site occupancy. This can be 
achieved, for instance, via a 2D Mott insulator state [22] . 
Subsequent patterning of the system can then provide 
the necessary spatial access to further manipulate and 
read-out the system by high-fidelity optical techniques. 

The depth of the optical lattice in Fig. [2] is 18 recoil 
energies (E r = tt 2 h 2 / (2ml 2 ) , with m being the rubidium 
mass and I = 600 nm). For these parameters the sin- 
gle particle tunneling time t = h/(SJ) amounts to 20 ms 
[23] . Here, h is Plancks constant and J ist the tunneling 
matrix element. As the preparation time (144 ms for the 
pattern in Fig.[2|c)) plus the subsequent imaging dura- 
tion (45 ms) is ten times larger than t, one might expect 
that a substantial fraction of the atoms rapidly refills the 
sites after the preparation and no structure is visible in 
the images. However, we only find a small population 
in the addressed lattice sites (12% in average). A some- 
what higher fraction of 25 % is found for the single defect 
in Fig.[2|a) - in accordance with what one might expect 
as atoms from four neighboring sites can tunnel into the 
defect. Isolated sites are stable as well. The atom num- 
ber in the single site shown in Fig.[2|b) only decays to 
55 % of its initial value and for the chain and the ring 
(Figs.gc) and (d)) we find even higher values of 80 % - 
90 % for the remaining population. These results suggest 
that tunneling in the patterned gas is highly supressed. 

In order to explain this observation we have to look 
at the system in more detail. The radial (axial) os- 
cillation frequency in each tube is uo r = 2tt x 12 kHz 
(uo a = 2tt x 170 Hz). With a temperature of T = 100 nK 
and a chemical potential of \i = h x 2 kHz the gas in each 
tube is one-dimensional. While the atoms are radially in 
the ground state they occupy more than ten vibrational 
modes in the axial direction. The central line density 
amounts to n\jj — 15 \ivn~ 1 and results in an interac- 
tion strength of 7 = 0.08. The dynamical evolution of 
the system after its preparation therefore corresponds to 
that of an array of coupled, weakly interacting ID Bose 
gases. An atom tunneling from a filled site with TV atoms 
into an empty site has to overcome the total interaction 
energy NU where U is the on-site interaction energy. 
As NU equals the chemical potential /i, such a tunnel- 
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FIG. 3: Single-shot image of the central part of the cloud. 
Each dot corresponds to a detected ion (800 in total). The 
integrated data along both lattice axes are shown on the side 
(black line). A fit to the data (red line) is used to determine 
the position of the lattice with an uncertainty of ±30 nm. 



ing process is highly off-resonant and cannot be compen- 
sated by the small tunneling coupling J. Consequently, 
the tunneling dynamics is suppressed. While tunneling 
of all N particles together is energetically allowed, it is 
highly unlikely as it requires a TVth-order tunneling pro- 
cess. In fact, the tunneling dynamics in the patterned 
gas has similarities with nonlinear self-trapping [24] and 
the observation of repulsively bound pairs [25]. However, 
as the axial level splitting is much smaller than the inter- 
action energy, tunneling into vibrationally excited axial 
modes of an empty tube might occur as an additional 
relaxation process. Moreover, as each tube can be de- 
scribed in local density approximation with an effective 
chemical potential fJL G s(z) = ji — \muj\z 1 , where z de- 
notes the direction along the tube, the above reasoning 
fails close to the edges of each tube. These effects might 
explain why we observe a fraction of atoms refilling the 
empty tubes. We plan to extend these studies in the fu- 
ture by changing the lattice geometry and preparing the 
tubes perpendicular to the electron beam. This allows to 
additionally measure the axial extension in each tube and 
to analyze thermalization processes during the tunneling 
dynamics. 

The optical lattice can exhibit short and long term 
drifts in both axes. We account for this by phase analyz- 
ing each individual image. In Fig.[3] we show a single-shot 
image together with the integrated linescans in both lat- 
tice directions. From a fit to the linescan we extract the 
actual phase of the lattice and displace the image accord- 
ingly prior to the summation. We routinely correct all 
images applying this procedure. In order to illustrate the 
drift motion of the lattice we show in Fig.[4^a) and (b) 
the phase distribution of a sequence of 300 images. While 
one axis is relatively stable over the whole sequence (cor- 
responding to a total measurement time of 1.5 hours) the 
second axis shows an almost uniform distribution. This 
becomes even more evident if we plot the phase evolution 
of each lattice axis as a connected chain of unity vec- 
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FIG. 4: Drifts of the optical lattice. (a,b) Appearance of lat- 
tice positions (translated into a phase between and 2tt) for 
300 experimental runs, corresponding to 1.5 hours of measure- 
ment time. The lattice vector is parallel (a) or perpendicular 
(b) to the condensate axis. (c,d) Illustration of the drift dy- 
namics for the data in (a) and (b). Each lattice position is 
converted into a unity vector whose direction is given by the 
corresponding phase. All 300 vectors are successively plot- 
ted. A perfectly stable lattice would result in a straight line. 
Small long term drifts appear as large circles, while short term 
fluctuations lead to a pronounced zigzag shape. 



tors, where the direction of each unity vector represents 
the actual phase of the optical lattice (Fig.Qc) and (d)). 
We attribute the short term fluctuations of the lattice to 
mechanical vibrations of the setup and air turbulances 
in both arms of each lattice axis. Long term drifts are 
probably due to thermalization processes of the setup. 
We have checked that drifts of the field of view due to 
electrostatic charging in the electron column play a minor 
role. The observation of larger fluctuations and drifts for 
one axis is compliant with a less stable opto-mechanical 
layout of this axis. 

For the addressing scheme, the situation is more in- 
volved as a motion of the optical lattice can lead to an 
error in the addressing protocol which cannot be com- 
pensated for in a post processing procedure. This is the 
case, for instance, if the electron beam points in between 
two sites, thus affecting both of them. Due to the high 
spatial resolution the addressing procedure even works 
for some error in the lattice position. For the images in 
Fig.[2]we have accepted all phases in an interval ± 0.6 x n. 
A phase of tt corresponds to the extreme case where the 
electron beam points in between two sites. With this 
acceptance interval, we have to discard on average 65 % 
of the images. Imperfect addressing also leads to longer 
depletion times. The dependence of the depletion time 



4 



0.8 



3 0.6 

(/) 

c 

o 

o 

2 0.4 



c 
o 

% 0.2 

Q. 

CD 
Q 

0.0 




d=90 nm 



0.0 0.2 0.4 0.6 0.8 1.0 

Addressing time (ms) 



50 100 150 200 250 300 
Distance d from site center (nm) 

FIG. 5: Depletion time constant for different positions of the 
electron beam with respect to the lattice site. At all distances 
we find an exponential decay. For central hits, the population 
decays with a time constant of 100 /is. The inset shows the 
depletion at a distance of 90 nm. The solid line is a fit with an 
exponential decay. The offset of 10 % represents atoms that 
are refilling the lattice site between the end of the preparation 
and the instant of imaging. 



constant on the position of the electron beam is shown in 
Fig.|5j Whereas for central hits we find a time constant 
of 100 /is, it is more than 10 times larger if the electron 
beam points at the edge of a site. For all positions we 
find an exponential decay of the atom number. In terms 
of writing time, a lattice site can be prepared in about 
one millisecond - provided that the site is hit in the cen- 
ter [26]. If we allow for a preparation time of 100 ms a 
total of 100 sites can be emptied. This is sufficient to 
tailor a large part of the gas. 

The ability to remove atoms from specific sites has a di- 
rect possible application for the study of quantum phases 
in optical lattices. In many cases, the temperature of the 
gas after loading in the lattice potential is too high to 
reach the ground state of the system. Therefore, further 
cooling in the lattice is necessary. One possible cooling 
strategy relies on the spatial separation of the gas into re- 
gions of high and low entropy which are subsequently de- 
coupled by a potential barrier from each other [27 j. Our 
technique would directly allow to remove atoms from re- 
gions of high entropy without the need to decouple them 
from the remaining system. 

In conclusion we have shown that arbitrary patterns of 
mesoscopic atomic ensembles in an optical lattice can be 
created and imaged by means of scanning electron mi- 
croscopy. The patterns are remarkably stable even for 
a moderate lattice depth and allow for in situ studies 
of tunneling processes in optical lattices. The technique 
provides large flexibility and can be applied to any atomic 
species, mixture of different species or molecules. The 
tailored quantum gases can be further processed by stan- 



dard manipulation and detection techniques and consti- 
tute a novel resource for cold atom physics. Ultimately, 
they might allow for the implementation of new schemes 
for quantum optical applications, quantum simulation 
and quantum information processing. 

This work was supported by the DFG within the 
Emmy Noether-program and the University of Mainz. 
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